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The self-consistent random-phase approximation (SCRPA) is reexamined within a multilevel- 
pairing model with double degeneracy. It is shown that the expressions for occupation numbers 
used in the original version of SCRPA violate the particle number for non-symmetric particle-hole 
{ph) spectra. A renormalization is introduced to restore the particle number, which leads to the 
expressions of occupation numbers similar to those derived by Hara et al. for the ph case. The results 
of calculations within the ph symmetric case show that this number-conserving SCRPA yields the 
energies of the ground state and first excited state of the system with n -I- 2 particles relative to 
the ground state of the system with particles in close agreement with those obtained within the 
g . original SCRPA. However it gives a slightly larger correlation energy. 
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The random-phase approximation (RPA) has been a powerful tool in the theoretical study of many-body systems 
such as atomic nuclei. An essential ingredient of the RPA is the use of the quasiboson approximation (QBA), which 
\ considers fermion pairs as boson operators, just neglecting the Pauli principle between them. Within the QBA a set 
^-H ■ of linear equations, which is usually called the RPA equation, is derived, which makes computationally demanding 
J> ' problems become tractable. However, because of the violation of Pauli principle within the QBA, the RPA equation 
, breaks down at a certain critical value of the interaction's parameter, where the RPA yields imaginary solutions. 

Several approaches were developed to remove this inconsistency of the RPA. One of the popular ones is the renor- 
2 ' mahzed RPA (RRPA) 0,111111. The RRPA includes in the expectation value over the ground state the contribution 
I of the diagonal elements of the commutator between two fermion-pair operators. In this way it takes the Pauli princi- 
^ pie into account approximately. This includes the so-called ground-state correlations beyond RPA, which eventually 
I renormalize the interaction in such a way that the collapse of RPA is avoided. However, the tests carried out within 
exactly solvable models also showed that there is still a large discrepancy between the solution obtained within the 
RRPA and the exact one beyond the RPA collapsing point (See Ref. 3 e.g.). 

The situation has been significantly improved recently within the self-consistent RPA (SCRPA) [10,[3 due to the 
O , inclusion of screening corrections in the SCRPA equation. These screening corrections are in fact the expectation 
values of the products of two fermion pairs in the correlated ground state. As the result the sign of the interaction is 
reversed so that, within a particle-hole {ph) symmetric multilevel-pairing model with double degeneracy (the so-called 
picket-fence model), the SCRPA yields the solutions very close to the exact ones for the correlation energy of the 
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system with 57 particles, as well as the energy of the first excited state of the system with + 2 particles 

Realistic nuclear single-particle spectra are in general ph non-symmetric, which means that the particle-particle 
' {pp) submatrix A and hole-hole {hh) submatrix C of the pp-RPA equation do not have the same dimension. The 
asymmetry is particularly strong, e.g., in light neutron-rich nuclei (g, for which the effect due to Pauli principle 
cannot be neglected. It is, therefore, worthwhile to reexamine carefully the SCRPA before applying it to realistic 
nuclei. 

Thepresent paper employs the same picket-fence model, which was used to tested the validity of the SCRPA in 
Refs. HQ- It will be shown that, in the general ph non-symmetric case, using its original expressions of ground-state 
correlation factors, the SCRPA violates the particle number. A simple and consistent way to restore the particle 
number will be introduced and the consequences will be discussed. 

The paper is organized as follows. The outline of the SCRPA for the picket-fence model is presented in Sec. 
The violation of particle number within the SCRPA for ph non-symmetric case and the construction of a number- 
conserving SCRPA are discussed in Sec. Illll The results of numerical calculations are analyzed in Sec. IIVI The paper 
is summarized in the last section, where conclusions are drawn. 
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II. SCRPA EQUATION FOR THE PICKET-FENCE MODEL 

The detail derivation of self-consistent pp RPA equation, which is simply called SCRPA equation hereafter has been 
described Refs. [ESQ- The present section recuperates only the brief outline of the SCRPA for the picket-fence 
model, which is needed for the discussion in this paper. 

A. Model Hamiltonian 

The picket-fence model consists of two-fold equidistant levels interacting via a pairing force with a constant 
parameter G. The model Hamiltonian is written as 

a n 
7J = ^(e,-A)7V,-G^ i^tp^- , (1) 

where the particle- number operator A^^ and pairing operators p/, Pi are given as 

= c\c, + cL,c_. , P} = c\cl^ , = (P;)t . (2) 

The exact commutation relations between the operators Ni, Pj, and Pi are 

[P,,p]]^5,j{l-Ni) , (3) 

[iV„ Pj] = 2%Pj, [N,, P,] = -2%P,- . (4) 

The single-particle energies take the values ti = ie with i running over all levels. The original version of the SCRPA 
in Refs. 0,13 applied only to a p/i-symmetric spectrum, where there are as many particles as levels (half filling). 
This means that, in the absence of interaction (G =0), the lowest fi^ = 17/2 levels are occupied with N = Vl particles 
(two particles in each level). However, in general, the number Q.h of hole levels is not necessary to be the same as the 
number Vtp of particle levels, i.e. Vlh ^ ^ ri/2. Numerating particle {p) and hole [h) levels from the levels closest 
to the Fermi surface, the particle and hole energies are equal to €p = e{D.h + p) and eh = ^{^h — ^ + 1), respectively, 
with h indices running from 1 to fi/i, and p indices running from 1 to f2p = 17 — 0,^. The Fermi level A is defined in 
the middle of the the first h and the first p levels, i.e. 

+ (5) 

Therefore, in the ph symmetric case [Vlp = Vlh = 17/2), the Fermi level A = [e{il + 1) — G]/2 is located in the middle 
of the single-particle spectrum. 
Using the notation of Ref. 

Mp^Np, Ah Ql = Pl , Qh = -PI , (6) 

and also introducing the ground-state correlation operators Dp and Dh 

Dp^l- Mp^l- Np , Dh = l- Ah , (7) 

the exact commutation relations lO and Q can be transformed into 

IQp^Qp'] = Spp'Dp , [Qh,Qh>] — ^hh'Dh , (8) 

[Ah Q]] = 2%Q] , [Ah Qj] = -26,jQ, . (9) 

Using the definition (jSJ of the Fermi energy A together with notations and Q, Hamiltonian Q is rewritten in 
the following form 

H=-enl + Y^[eip^^) + ^]Ah + Y.Hh-l) + ^]AIn 

p=l h=l 
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(10) 



hh' 



ph 



which, in the ph symmetric case, coincides with Eq. (13) of Ref. 



B. SCRPA equation 

The SCRPA equation is derived based on the RRPA additional and removal operators, which have the form 



4 = E^p"^?i-EW^ 



(11) 



and 



(12) 



respectively , with the abbreviation 



o 



i ^ p, h 



(13) 



to denote the renormalized operator of an operator Oj . Operator Aj^ transfers the states in a system with il particles 

to those of a system with f2+2 particles. Operator R\ transfers the states of an fi-particle system to those of a system 
with particles. The brackets (...) = (fi, 0| . . . \ ft, 0) denote the average over the correlated ground state 0) of 
the system with D, particles, which is defined as the vacuum of operators and R\, i.e. 



A^\n,o) = Rx\n,o) = o . 



(14) 



Using the exact commutation relation |(SJ) and the definition 114|l of the RPA ground state 0), one can see that the 
additional and removal operators satisfy the boson commutation relations in the ground state \Q,0) 



if the amplitudes X and Y satisfy the following normalization (orthogonality) conditions 

E ^p^p' - E ^kY/:' = 5,,' - E ^^^'^ - E ^i^Yp - -^^A' , E ^p^p' - E ^^^li = o 



p h 

while the closure relations 



E ^P^P' - E YpYp' - ^pp' ' E ^hXt^ - E ^/rn^ = ^h,, , E ^^^p - E ^p = » 



(15) 
(16) 

(17) 



guarantee the following inverse transformation of Eqs. Ullfl and (|12|1 



Y^x>:rx + Y.^I:ai 



(18) 



The SCRPA equation is obtained in a standard way by linearizing the equation of motion. The matrix form of the 
SCRPA equation for the additional mode is 



A B 
~B C 



yP 



yp 



(19) 
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where the submatrices A, B, and C were derived in Ref. using the definition Q as well as the exact commutation 
relations © and Q as 



(20) 



i?,, = -([ Q„ [ff,Q,]]) = G ^^^y , (21) 



(22) 



The expectation values of the products of two pair operators at the right-hand side (rhs) of Eqs. (|20|l and (|22|l are 



{QlQp' ) = {P^Pp' ) = J{D.p){Dp, ) Yp%) , (23) 



A 



{QpQh) = {QIqD = -{PlPp) = -{PlPh) = ^{Dp){Du) X>lYp^ , (24) 

A 

{QlQh') = V{Dh){Du')Y.Y/:Y,^ - {Pl,Ph) - 5h,ADu) , (25) 

M 

With {Pl.Pj,) = ^{Du){Dh-)Yx>:xii, . 

A 

They were derived using the inverse transformation (|18|l and the definition of the ground state (|14|l . The RRPA 
equation was obtained from Eqs. p9|) - (|22|l by using the factorization 

{D,D,) c {D,){D,) , (26) 
and neglecting all the expectation values {Qp,Qp), {QpQh), and {Q^Qh')- The RRPA submatrices have then the form 

Mp-\) + %]^pp'~G^iDp){Dp,) , (27) 



RRPA 



4HK 

pp 



nRRPA 
^ph 



G^{Dp){Du) , (28) 



1. 



C^^''^ = -2 [e{h - 2 ) + ^^^h' + GV{Dh){Dh') • (29) 
The RPA submatrices are obtained from the RRPA ones by putting Dp = Dh = 1, namely 

A^PA = 2[e(p-i) + ^]V-G, (30) 
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B^h^ = G , (31) 



Cf^,^^-2[e{h-\) + ^]5u,+G . (32) 



By using definition @, Eqs. (|23|l - l|25|l . and recalling that 



e^-A=[e(p-i) + ^] , (33) 



e.-A = -6(/.-i) + | = -[e(/i-i) + |] + G, (34) 



one can rewrite Eqs. H2()|l - H22|l in the notations of Ref. as ^ 

A„, = ([Pp, [H,p\A]) = 2V(V - A) 



+2GV^--"<^--''^-^+fi"^^-^-^ - ^ ((l^Mi=^ , (35) 

5,. = ([P„ [H,P{]]) . G , (36) 

V(l-(A^p))((iV,)-l) 

Cw,' = -([n, - 2(5w,'(e/,' - A) 



V((^'>) - l)((^v> - 1) 
C. Correlation energy 

The correlation energy E^^^^^ is defined as the difference between the energy £^ = {H) in the ground state |fi, 0) 
(|14|1 and the Hartree-Fock (HF) energy £hf- The former is easily obtained from Eq. (|10|l while the later is — ei^^- 
The final expression for the correlation energy is obtained as 

i?™ ^ {H) - £hf = E [e{p - ^) + f ] (1 - {D,)) + E Hh - i) + f ] (1 - {D,}) 

p h 

-G[Y,{QlQp') -'^T.^^l^i)^ ' (38) 

pp' hh' ph 

where {Q\Qi) and (Q^Qj,) are given by Eqs. ^ - and Di by Eq. g7|). 



There are several misprints in Eqs. (11) of Ref. namely the factor 2 in front of all (NiNj) in the numerators of the last terms at the 
rhs of submatrices A, B and C should be eliminated. Also, the factor 1 -- {Npi ) in the denominator of the last term of Bp/^ should be 
replaced with (Nf^) — 1, and the sign "— " in front of "i&hh' {^hh' ~ ^) in the expression of C^^' should be reversed. 
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For comparison, the correlation energy within the RRPA is derived here by approximating the Hamiltonian as 
H ~ i/i^i^PA = £^'^ (RRPA) + ^ El'^^'^AlA^ + ^ E^''''^r\r^ , (39) 

where £^^^ (RRPA) is the energy in the RRPA ground state, while the eigenvalues E™^^ and E™'^^ are the excitation 
energies (|42|l and H43|l of the additional and removal modes, respectively, which are obtained by solving the RRPA 
equation, i.e. Eq. (|19|l with submatrices (|27|l - H29|) . The correlation energy 

i;RRPA = £:f2(RRPA) - ^HF is obtained 
by calculating the expectation value of Hamiltonian (|39() in the HF ground state |HF). By using Eqs. and 112() 
as well as the definition of |HF), for which Pp|HF) = (HF|Ppt = P^^|HF) = (HF|P/, = 0, we finally obtain 

Ef]^r^ = fo"(RRPA) - (HF|iJ|HF) ^ E ^ " E E ^ ' (^0) 

fJ. h X P ^ 

The RPA correlation energy E^^.^ is recovered from Eq. (|40|l putting D,; = 1, namely 

Ef^^ ^ Yl E(^'n^ - E E^^EiYpf ' (41) 

II h A p 

with E^^^ and i?^^^ denoting the RPA excitation energies of the additional and removal modes, respectively. 



III. PARTICLE-NUMBER WITHIN SCRPA 



The SCRPA Eqs. H19|l - ()22|1 has flp solutions for the additional mode with positive eigenvalues 

^M=^M^'-4'>0, M=l,---,f^p , (42) 

which are excitation energies of the + 2 system relative to the ground state of the O system. Since the SCRPA 
equation for the removal mode has exactly the same form as that of Eqs. H19|l - H22|l with the only difference that 
{X^,Y/^, E^) should be replaced with {Yp,X^, —E\), the flh negative eigenvalues 

Ex = £'t^ - 4' < , A = 1, . . . , O,, , (43) 

of Eqs. (|19|l - (|22|l have physical meaning as excitation energies of the 51 — 2 system relative to the ground state of 
fl system. This set of equations should be solved self-consistently with the normalization conditions IjlGI) and the 
equations for the factor {Dp) and (-D/i), which represent the ground-state correlations beyond the RPA. It is clear from 
Eq. (Q) that, in the absence of ground-state correlations (beyond RPA), the ground state |O,0) becomes the RPA 
ground state |RPA}, for which (RPA|A|RPA} = {RF\D,\RF) = 1 (i = p, h) due to the QBA, where HF denotes the 
HF ground state. This means that (RPA|Afp|RPA) = and (RPA| A^/JRPA) = 2. In this case, the SCRPA equation 
reaches its RPA limit with the RPA submatrices H30|l - H32|l. In the general case, < (Di) < 1 {i = p, h) since 
< {Np) < 1 and 1 < {Nh) < 2. 

A. Violation of particle number within SCRPA 

In order to derive the equations for the factors (Di) (i — p, h) Refs. 0,0 employed a procedure similar to the one 
used in Ref. with the representation 

= 2PIp, , (44) 

which becomes exact for the picket-fence model. Using Eqs. (Ujl, H23|l and (|25|l . one finds immediately from the Eq. 
(EH) 

(Np) = 2{Dp) J2iY^'r = 2(1 - {Np}) J2iY^')' , 

A A 



(Ar,)=2(D,)[l + ^(y,n'] = 2((A^/.)-l)[l + E(^^'')'] • (45) 



7 



This yields 



{N,) = 1-{D,) = 2{D,)Y,{Y^')\ 



{N,,) = l + {D,)=2[l-{D,)J2iY>:f] , 



(46) 



with 



{Dp) = 



1 + 2Ea(^p^)' 



(Df,) = 



1 + 2E,(W 



(47) 



This resuhs is a special (degenerated) case of the equations for the pp and hh ground-state correlation factors {Dppi ) 
and {Dfih') in the general realistic spherical shell-model basis, which is derived here using the general expression of 
the relation (|44|l in the form 



JMj' 



with 



4, ( JM) = ^ (jm/mV'M')4„.c],,„, , P,,. ( JM) = [4, ( JM)]t 

mm' 

Inserting in the rhs of Eq. 148|) the general expression for pL(JM) 



(48) 



(49) 



and using Eqs. (|14|l and (|15|) . the final equations for {Djji) = 1 — rij — n^/ is obtained in the form 



(^pp') = i-E(2>/+i)E 



J, A 



{Dpp") 



lyJA 12 

(2jp + 1) 



(^P"P') 



yJA 12 



(2jp' + 1) 



(50) 



h" 



lyJf^ 12 12 



(2j/. + 1) 



i2jk' + 1) 



(51) 



where jp (jh) denotes a p (h) orbital angular momentum, and J is the total angular momentum (multipolarity of the 
excitation) . Obviously, Eq. (|T7jl is recovered from Eq. (|5T|l in the degenerated case, when J = jp = jh = , p = p' , 
and h = h' . 

Equations (gHl) and are the result given by Eq. (13) of Ref. 6] 2, and Eq. (30) of Ref. This result for the 
pp case is similar to what obtained previously in Ref. 4], but for the ph case, according to which 



{{Np))^2j2{{Dph))iY;,f , {{Nh))^2[l-J2{{Dpk))iY;Hf] , 



(52) 



^ The index A in the sum at the rhs of the expression for {Np) in Eq. (13) of Ref. Q has been misprinted as /i, although this did not 
affect the results of calculations for the ph symmetric case, for which fi = X. 
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where 

Here, to avoid confusion with the notation for the pp case, the double brackets ((. . .)) are used to denote the av- 
erage over the correlated ground state with respect to the ph renormalized RPA operators. Except for this formal 
similarity, the essential difference between the ph and pp cases is that Eq. (|52f) for the ph case always conserves 
the particle number, while Eq. for the pp case, in general, does not. Indeed, in the ph case, Eq. gives 

J2p{{^p))J2p~^J2hii^h)) — '^Y^h'^ 1 = because the two sums at the rhs of Eq. cancel each other exactly 

and, therefore, independently on how {{Dph}) is estimated. Meanwhile, in the pp case, Eq. in general, violates 

the particle number because it gives 



= + 2^ [Y,{D,){X^)^ - Y.{D,MY/,r] -2Y^{D,) ^ n , (54) 

MP h p 

as 

Y,[Y,{D,){Xp'-Y,{DH)iY,^n^Y.(^p) (55) 

p. p h p 

unless the condition ~ \Y/i^\ is assumed, which means {Dp) = {D^). This condition is satisfied only when the 
full ph symmetry holds, i.e. Op = 0/j = 0/2 for the equidistant spectrum. In the general ph non-symmetric case, 
i.e. when Op ^ O/j ^ 0/2 and/or the spectrum is not equidistant, there is no normalization condition such that H55|) 
becomes an equality since this would be incompatible with the normalization condition (|16|l for the SCRPA X!^ and 
amplitudes. 

The measure 50 of particle-number violation can be estimated by expanding {Di) into the power series of {Y^)"^. 
By using the normalization and closure relations H16() and (|17|l . the lowest order of this expansion yields 

m ^ I Y.{Np) + Y.{Nu) - 0| = 12^] [Y,{Dp){X!;f Y.^Du){Y>:f] ~2Y,{Dp)\ 



-|2E{E[i-2E(^p')']W-E[i-2E(^/r'n(^/r)'}-2^[i-2E(^p')']| (se) 

p p \ hp' p X 

= 4| Y.(Yn'(y,f)' - E(^p')'(^p'')'l - ^(^') > 

pp'h XX' p 

which means that the particle-number violation is expected to be small at least within the validity region of RPA, 
where \Y^'^\ are small. 



B. Restoration of particle- number conservation within SCRPA 

Equations (|46|) have been derived making use of Eq. 1)44(1 . which is compatible only with the exact ground state 
|O,0) (I14|) . However, as has been discussed in details in Ref. 0, such exact ground state does not exist within 
the SCRPA, except for the case with = 2, where the SCRPA and exact solutions coincide. Consequently, the 
SCRPA formalism still contains some violation of Pauli principle, which leads to the particle-number violation in the 
ph non-symmetric case. 

In order to restore the particle number within the SCRPA, let us notice that the essential point of SCRPA and 
RRPA is the renormalization of the operators and Qi ((HJ) in such a way that the renormalized operators incorporate 
the effect of Pauli principle due to their fermion structure, but behave at the same time like ideal boson operators with 
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respect to the expectation value (...) = 0| ... 0) The result of such renormalization yields the operators 

Qj = Q\I \J (Di) and Qj = Qij \J (-D,;} in Eqs. and p2|l . which satisfy the following exact relation 

([Q„Oj]) = *y . (57) 

This relation means that, as far as the calculation of expectation values is concerned, the commutator [Qj,(5j] can 
be safely replaced with its ground-state expectation value ([Qj,Q]]), namely 

n.M\\^m.M^)^k, . (58) 

i.e. and Qj are now considered as ideal boson operators (without fermion structure). From now on the derivation 
is proceeded only with expectation values using the replacement (|58ll under the condition 

(4A,,.) = (i?li?v)-0 (59) 

instead of the vacuum condition H14|l.^ Using Eq. 0, one can see that the renormalized operators and satisfy 
the following exact commutation relations with operators Mi : 

\M,,q\\ = 2J,,q1 , [M„Qj] = -2,5,,Q, . (60) 

Since the standard derivation of RRPA equations is based on the algebra (|58|l and (|60|l in terms of the boson operators 

and 0,0, 0|j in order to derive the equations for the renormalization factor (Z?;), the fermion operators Mj are 
also bosonized so that Eq. (|60|l still remains intact. Such boson representation exists and equal to 



M, = 2q1q, , (61) 



which fulfills Eq. H60() exactly since 



[M„Q]] = 2gI[g„Q]] = 2%gl , 

[M„Q^.] - 2[qI,Q^.]Q^ - -2%Q, , (62) 
due to Eq. H58|) . Representation ()61|l is apparently different from Eq. (|44|l since the latter is equivalent to 

M, = 2Q|g, = 2(A)qIq, (63) 

due to definition The commutation relations between Eq. (|63|l and operators and are different from the 
exact relations (|60|l by the factor which causes the particle-number violation discussed in the preceding section. 
Using the representation (|61|l instead of ^S\\\ together with Eqs. (|23|l and H25(l immediately leads to 

(7V^) = (M,) = 2^(y/)2 , (iV,) = 2 - (M,) = 2[1 - Y.^Yi: f\ ' (^4) 

A 

or 

{D,) = 1 - 2^(r/)2 , (Z?,) = 1 -2Y^iJi:f , (65) 

A M 

instead of Eqs. (|46() and (|47|1 . As has been mentioned previously, since 0< (A'p) < 1 and 1 < (A^h) < 2, the values of 
ground-state factors given in Eq. (65) should also satisfy < < 1. These results are similar to what obtained by 



^ This situation is somewhat similar to that of the statistical formalism, where a quantum mechanical ground state |0{/3)) {/3 is the inverse 
temperature) so that {0(/3)|C'|0(/3)) = Tr{C'I'} {T> is the density operator) does not exist. The expectation value (O) is then replaced 
with the statistical average over the grand canonical ensemble with D = e~'^(^~^^)/Tr{e~''(^~^^) }. 
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FIG. 1; Particle-number violations within SCRPA for ph non-symmetric cases with f2 = 10 and ^Ih ~ 1 (solid line), 2 (dashed 
line), 3 (dotted line), and 4 (dash-dotted line) as functions of interaction parameter G (in units of the level distance e). 

Hara in Ref. and, later by Rowe in Ref. but for the ph case. Therefore, the version of SCRPA (RRPA), which 
uses Eqs. (jMIl and (jHSJ instead of Eqs. and (gZIl, will be referred to as Hara SCRPA (Hara RRPA) hereafter 
That Eq. I|64(l conserves the particle number is straightforward, making use of the normalization and closure relations 
(|16ll and (|17|) . Indeed, using Eq. Ht)4|l instead of Eq. H4t)|) . one finds in the same way as it was done in proving the 
particle-number conservation within the pp RPA that 

p h P A h 11 

= n+2{ E[E(^p")' - 1] - E E(^.r)'} (66) 

P fJ. h fJ. 

= ^+Y. [E(^p )' - E^^/T)'] -% = ^ + %-%^^. 

MP h 

Therefore, the Hara SCRPA and Hara RRPA always conserve the particle number exactly. 

IV. NUMERICAL ANALYSIS 

The calculations were carried out for several values of and e = 1 MeV within SCRPA and RRPA. The most 
representative case with f2 = 10 is selected here for discussion. For simplicity, the factorization (|26|l was used, which 
has been verified in Ref. J] to yield excellent results compared with those obtained when an involved set of nonlinear 
equations for the expectation values (MiMj) was solved instead. This factorization does not affect the discussion 
regarding the particle-number restoration in this work. 

A. Degree of particle-number violation within SCRPA 

Shown in Fig. ^is the quantity SQ,/{2ilh) as a function of the interaction parameter G (in units of level distance e), 
which has been obtained within SCRPA for ph non-symmetric cases with the number of hole levels f2/i — 1, 2, 3, and 4. 
The particle-number violation increases with G and with the asymmetry of the ph single-particle space. The strongest 
violation of about 2% is observed at the strongest asymmetry, i.e. with il/j = 1 and Q.p — 9 (solid lines), at the largest 



* The expressions for the factor Dppi and Df^f^i within the Hara-SCRPA for the general shell-model spherical basis are obtained from Eqs. 
1511 putting Djji at the rhs of Eqs. 1511 equal to 1. 
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value of G = 0.45 MeV shown in the figure. In all other cases plotted on this figure, the particle-number violation is 
smaller than 1%. With increasing fJ^, the symmetry is gradually restored, and the particle-number violation decreases 
to reach zero at fi/j = 5. Results of our calculations for larger fl also show that the particle- number violation within 
SCRPA decreases with increasing the particle number. 



Shown in Fig. [21 arc the correlations energies of the system with — 10 particles, as well as the energies Ei and 
E2 of the ground state and first excited state, respectively, of the system with ft + 2 = 12 particles relative to the 
ground state of the fJ-particle system as functions of the interaction parameter G (in units of e) . They were obtained 
within the RPA, RRPA, SCRPA, and are plotted in comparison with the exact results. The RRPA gives a quite 
good description of the correlation energy, which practically coincides with that given by the SCRPA and the exact 
result for G < 0.45 MeV. However, the RRPA fails badly in describing the the ground-state Ei and first-excited-state 
E2 energies of the -I- 2 system. Here, although the RRPA results do not collapse at Gcr — 0.34 MeV as the RPA 
results do, they decrease monotonously, while the exact results as well as those given by the SCRPA increase with 
increasing G (Cf. Refs. The results obtained within the Hara-SCRPA are close to those given by the RRPA, 

but fail to converge in this model at G > Gcr- The Hara-SCRRPA, which conserves the particle-number exactly in 
ph non-symmetric cases, offers very close results to those given by the SCRPA for the Ei and E2 energies within 
the whole interval of values for G under consideration. However, the correlation energy i?corr obtained within this 
number-conserving version of SCRPA is slightly larger than the exact result, and the discrepancy is clearly visible 
already starting from G > 0.3 MeV. 

In general, the feature depicted in Fig. |21is similar to that of the ph case considered in Ref. 4], where the solution 
obtained by using Eq. H53|l approaches the exact solution at large G, while the one offered by the Hara approach 
fails to describe it, never approaching zero. This result comes from the overestimation of ground-state correlations 
beyond RPA within the Hara approach, which can be clearly seen by examining the ground-state correlation factors 
(Di) and/or the occupation number Ni. The factor (Dp), which is the same as (Dh) for the symmetric case of the 
picket-fence model, is shown in Fig. 13 This factor decreases from 1 with increasing G, approaching zero as G ^ 00. 
The deviation from 1 is stronger at the level closer to the Fermi one. The difference between the results obtained 
by using Eqs 1)53(1 (the SCRPA) and 1(6 5|l (the Hara-SCRPA) is strongest for the lowest particle level, in which the 
Hara-SCRPA gives stronger ground-state correlations beyond RPA. It also increases with increasing G in line with 
the results obtained for the ph case in Ref. Q . This also explains the larger discrepancy between the two approaches 
in the description of correlation energy £^corr, while the differences in the energies Ei of the ground states, and E2 of 
the first excited states are relatively smaller. 

The behavior of {Di) leads to the change of the occupation number {Np) and {Nh) as shown in Fig. 0] As small 
G, the function Ni approaches the stair case with (Np) ~ and (Nh) — 2 for all ei. As G increases, the deviation 
from the stair case becomes more and more evident with the decrease of (Nh) from 2 and increase of (Np) from 0. At 
G/e = 0.4, e.g, (Nh) becomes 1.85 and (Np) reaches 0.4 for levels closest to the Fermi one. The deviation caused by 
the Hara-SCRPA is always stronger than that given by the SCRPA. Since (Di) — * as G ^ 00 we have from Eq. I(t)5() 
(the Hara-SCRPA) the sum Ea(^p^)^ = J2^,(Xh? 1/2, which make Np ^ Nh ^ 1. For the SCRPA at the value 
J2xiYp^f = llp,{Yh? = 1/2 one obtains {Dp) = {Dh) = 1/2, which lead to {Np) = 1/2 and {Nh) = 3/2. Again, this 
shows that ground-state correlations beyond RPA are stronger within Hara-SCRPA than within SCRPA. 

The exaggeration of the ground-state correlations beyond RPA within the rcnormalization procedure, which leads 
to the number-conserving (Hara) type expressions I(ti5() was pointed out before by Rowe in Ref. 2], where, by using the 
number-operator method to insert the number operator twice at the center of (A^i), he found that the ph ground-state 
correlation factor {{Dph)) became J^i, l^phl'^ instead of 2J2u l^phl'^- The result of an infinite expansion by inserting 
repeatedly the number operator at the center of {Ni) is not available for ppKPA at this stage. However, the observation 
by Rowe suggested that the real {Dp) and {Dh) might be closer to 1 than those given by Eqs. l(^ . Therefore, a test 
was also carried out here by parametrizing {Dp) and {Dh) within the Hara-SCRPA to see if it is possible to achieve 
results as good as those given by SCRPA for all three quantities Econ, Ei and £'2. For this test, we used 



instead of Eqs. ((65(1 and repeated the calculations. The result of this test shows that the values Econ, Ei, and E2 
obtained within the SCRPA can be fitted simultaneously rather well within such "parametrized" Hara-SCRPA with 
the parameter a ~ 1.9 . These results are shown in Table^in comparison with the SCRPA ones. Such "parametrized" 
Hara-SCRPA also conserves exactly the particle number in the ph non-symmetric case as the Hara-SCRPA does. 



B. Correlation, ground-state and excited-state energies 




a < 2 , 



(67) 
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FIG. 2; Correlation energy Ecoir in the system with SI = 10 particles (a), ground-state energy Ei (b) and first-excited-state 
energy E2 (c) of the system with f2-|-2 = 12 particles relative to the ground-state of Jl-particle system as functions of interaction 
parameter G obtained within RPA (dotted line), RRPA (dash-dotted line), Hara-RRPA (thin dashed line), SCRPA (thick solid 
line), Hara-SCRPA (thick dashed line), and exact (thin solid line) calculations. 



V. CONCLUSIONS 



The present work shows that the SCRPA violates the particle number in the ph non-symmetric case if the occupation 
numbers are calculated according to Eq. H46|l for the picket-fence model, which is a limit of Eq. (|51|l for the general 
shell-model spherical basis. Within the ph non-symmetric picket-fence model this particle-number violation increases 
with the asymmetry and interaction strength G, but it decreases with increasing the particle number. However, within 
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FIG. 3: Ground-state correlation factor (Dp) for particle levels as functions of interaction parameter G. The solid and dashed 
lines denote results obtained by using Eqs. H47|l and 1)65^ . respectively. The numbers next to the lines numerate the particle 
level, to which the lines correspond. 
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FIG. 4: Occupation numbers (Nh) for hole levels (a) and (Np) for particle levels as functions of single-particle energies at 
several values of G/e indicated as the numbers next to the circles. The full and open circles denote results obtained by using 
Eqs. 1471 and 165L respectively. The lines are drawn to guide the eyes. 



TABLE I: Correlation energy Ecow, ground-state energy Ei, and first-excited-state energy E2 obtained within the 
"parametrized" Hara-SCRPA using Eq. JSTJ with a = 1.9 (a) in comparison with those given by the SCRPA (b) for — 10. 
All the values are given in MeV. 







£1 


£2 


G 


a 




b 


a 


b 


a 


b 


0.01 


-0.3265 X 10" 


'i 


-0.3265 X 10~^ 


1.0000 


1.0000 


3.0000 


3.0000 


0.05 


-0.9018 X 10" 


2 


-0.8563 X 10~2 


1.0005 


1.0005 


3.0001 


3.0001 


0.10 


-0.3850 X 10" 


1 


-0.3652 X 10"^ 


1.0032 


1.0033 


3.0063 


3.0063 


0.15 


-0.9263 X 10" 


1 


-0.8808 X 10"^ 


1.0111 


1.0112 


3.0194 


3.0194 


0.20 


-0.1760 




-0.1686 


1.0278 


1.0279 


3.0811 


3.0460 


0.25 


-0.2930 




-0.2846 


1.0566 


1.0563 


3.0922 


3.0926 


0.30 


-0.4470 




-0.4424 


1.0985 


1.0970 


3.1645 


3.1657 


0.35 


-0.6401 




-0.6475 


1.1514 


1.1481 


3.2670 


3.2703 


0.40 


-0.8823 




-0.9023 


1.2118 


1.2058 


3.4011 


3.4087 


0.45 


-1.1608 




-1.2094 


1.2755 


1.2667 


3.5654 


3.5804 
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the interval of values for G under consideration (G < 0.5 MeV), we also found that the particle-number violation 
reaches at most around 0.2% for the most asymmetric case with the level number — 10, where the number of hole 
levels ilfi — ^, and number of particle levels — 9 . In all other less asymmetric cases this violation is smaller than 
0.1%. 

In order to maintain the exact particle-number conservation within the SCRPA, a rcnormalization was proposed, 
which represents the number operator in terms of the product of renormalized pairing operators. As a result, a 
number-conserving SCRPA was derived, which is called Hara-SCRPA as it has the equations for the occupation 
numbers similar to what obtained in the pioneering works by Hara et al. but for ph case [ij. The results of numerical 
calculations show that the Hara-SCRPA yields the ground state energy and energy of the first excited state of the 

+ 2 system very close to the corresponding values obtained within the SCRPA. However, the correlation energy, 
which the Hara-SCRPA offers, is slightly larger than that obtained within SCRPA. 

The results of the present study also indicate that, in realistic calculations using non-symmetric single-particle 
spectra within RPA, in particular for light systems, one should carefully examine the violation of Pauli principle to 
see if it is important to include the ground-state correlations beyond RPA. As a matter of fact, thepreliminary results 
of RRPA calculations, which were carried out recently for ^^'^"'Be using the Gogny interaction 9], have shown that 
ground-state correlations beyond RPA increased the correlation energy by 20 - 24% compared to the RPA results. 
This shifted up the ground state energy by 13% for ^^Be and 48% for ^^Be. At the same time the particle-number 
violation within the RRPA due to the use of Eq. H51|) did not exceed 0.2%. In this case the SCRPA can be still 
well justified, and has the advantage over the Hara-SCRPA as the former offers a better description of the correlation 
energy. 

In the cases where the particle-number violation cannot be neglected (e.g. > 1%) in calculations with realistic 
spectra and interactions a number-conserving approach like the Hara-SCRPA proposed in the present work might 
have to be used instead of the SCRPA. However, the improvement of the correlation energy in this case cannot be 
achieved by simply renormalizing RPA as has been done in the approaches under discussion. The test of parametrizing 
the Hara-SCRPA to yield all the three energies i?corr, Ei, and E2 close to the values given by the SCRPA suggests that 
higher-order correlations may have to be included in order to reproduce all these three quantities within a number- 
conserving SCRPA. This indicates coupling to configurations more complicated than the ph, pp, or hh ones should 
also be taken into account. 
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